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We study the dissipation physics of one dimensional mesoscopic superconducting quantum 
interference device array by using the field-theoretical renormalization group method. We observe 
length scale dependent superconductor-insulator quantum phase transition at very low temperature 
and also observe the dual behaviour of the system for the higher and lower values of magnetic field. 
At a critical magnetic field, we also observe a critical behaviour where the resistance is independent 
of length. 
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PACS numbers: 74.78.Na Mesoscopic and nanoscale system, 74.90-1- Other topic in superconduc- 
tivity, 05.10.Cc Renormalization Group methods 
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1. Introduction: The transition at zero temperature 
or at very low temperature belongs to the class of 
phase transition that is driven by the quantum fluc- 
tuatitons of the system The quantum fluctuations 
are controlled by the parameters of the systems such 
as the charging energies and Josephson couplings of 
the Josephson junctions array 2]. Here we present a 
field-theoretical renormalization group study to find the 
quantum dissipative phases of lumped superconducting 
quantum interference device (Squid). The quantum 
fiuctuations of the system are controlled by the ex- 
ternally applied magnetic field and a (= a ratio of 
quantum resistance to tunnel junction resistance). In 
the Squid total current of device is modulated through 
the applied magnetic flux. Therefore the total current 
in the Squid is / = 2IcSin{9)\cos{^)\, where Ic is the 

critical current, $ is the magnetic flux and $o(= 55) 
is the flux quantumi, 9 is the phase of superconducting 
order parameter. Similarly the Josephson coupling is 
also changing Ej = 2Ejqcos{9)\cos{^)\ , where Ejq 
is the bare Josephson coupling. So one can think that 
a lumped Squid system (Fig. lA) can be described in 
terms of an array of superconducting quantum dots 
(SQD) but with a modulated Josephson coupling (Fig. 
IB) and critical current [3,, JJ. This effective mapping 
will help us to analyze (analytically) the system in 
detail. The experimentalist of Ref.3 and Ref.4 have also 
considered the mesoscopic Squid system as a modulated 
Josephson junctions array. We have been motivated 
from these well-accepted experimental findings [3, 
The plan of this manuscript is the following. Section 
(II) contains the analytical derivations and the physical 
explanation for the occurrence of quantum dissipative 
phase in the lumped Squid system. Conclusions are 
presented in section (III). 



2. RENORMALIZATION GROUP STUDY FOR 
QUANTUM DISSIPATION PHASE IN THE 
MESOSCOPIC LUMPED SQUIDS 

We now present the basic dissipation physics of meso- 
scopic Squid systems. The source of dissipation at very 
low temperature is due to the appearance of phase slip 
centers. Phase slip centers are of two kind: one is ther- 
mally activated phase slip center valid near to the su- 
perconducting transition temperature, well described by 
LAMH theory 0, i, 0], and the other, quantum phase 
slip (QPS ) center occurs at T=0 K or at very low tem- 
perature due to the quantum mechanical tunneling in dif- 
ferent metastable states |8|. The most important type of 
fluctuation which occurs during this QPS process is that 
phase of the superconducting order parameter is chang- 
ing by ±27r at a point in the system and the amplitude of 
the order parameter vanishes at one point. Appearence 
of QPS in different low dimensional superconducting sys- 
tems is a common phenomena, we have cited only those 
references which are related to our problem [1, H, S S] ■ 
At flrst we prove the appearence of thermally activated 
phase slip centres in lumped Squid system. We consider 
a Squid and use the Ginzburg-Landau theory to illustrate 
the physics of QPS. A mesoscopic Squid is supercon- 
ducting ring (nano scale size) with two tunnel junctions. 
We therefore use the cylindrical polar coordinates. The 
superconducting order parameter inside the mesoscopic 
Squid is "0(0) = "00 e™"^, where (/) is the azimuthal angle, 
n is integer (winding number) and -00 is constant. Here 
the Squid is pierced by the magnetic flux $. The vector 
potential along the tangential direction is = 
{ R is the radius of the Squid). The free energy corre- 
sponding to this wave-function and vector potential is 
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First term is the free energy contribution in absence of 
magnetic flux, m* is the effective mass of the quantum 
system and /xg is the permiabiUty, K, is volume of the 
Squid, ci is a constant and the last term represents the 
vacuum magnetic field energy. Free energy is minimum 
when $ = ri$o- Depending on the externally applied 
magnetic flux, system is in one of the metastable min- 
ima and a persistent current flows around the mesoscopic 
Squid to maintain the superconducting state. The sys- 
tem can jump from one metastable minimum to the other 
minimum to lower the energy. This is the source of dis- 
sipation and the decay of persistent current. Such an 
event corresponds to a change in the winding number (n) , 
hence it is a phase slip. The rate of thermally activated 



phase slip is — 



''b'^ J where Cb ^ Vs 



is the 



barrier energy. For the bulk Squid, K, is large and hence 
the appearance of phase slip centers is quite unusual. But 
for the mesoscopic system of nano-scale size such phase 
slip centers are likely possible. In the nano-wire, there 
are many evidences for the presence of thermally acti- 
vated phase slip centers near the superconducting tran- 
sition temperature Q . The basic transport properties of 
low dimensional tunnel junction systems and nano-wire 
have many similarities [lo| . So we also expect the 
appearance of thermally activated phase slip centers for 
mesoscopic lumped Squid system. Temperature region 
close to the transition temperature of lumped Squid sys- 
tem 0, 3] is the region for thermally activated phase slip 
centers. 

Here we prove the appearance of QPS in SQD array with 
modulated Josephson coupling through an analysis of a 
minimal model. We consider two SQDs are separated by 
a Josephson junction. These two SQDs are any arbitrary 
SQDs of the array. Appearance of QPS is an intrinsic 
phenomenon (at any junction at any instant) of the sys- 
tem. The Hamiltonian of the system is 



H 
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where rii and Oi are respectively the Cooper pair den- 
sity and the superconducting phase of the i'th dot and 
C is the capacitance of the junction. The first term of 
the Hamiltonian present the Coulomb charging energies 
between the dots and the second term is nothing but 
the Josephson phase only term with modulated coupling, 
due to the presence of magnetic flux. We will see that 
this model is sufficient to capture the appearence of QPS 
in SQD systems. Hence this Hamiltonian has sufficient 
merit to capture the low temperature dissipation physics 
of SQD array. In the continuum limit, the partition func- 
tion of the system is given by, Z = ^ D9{x, r) Sq(x,t}^ 
where Sq ^ J dr J dx^ [{dr0{x,T)f + {d^e{x,Tyf]. 
The action is quadratic in scalar- field 9{x,t), where 
9{x, t) is a steady and differentiable field, so one may 
think that no phase transition can occur for this case. 



This situation changes drastically in the presence of 
topological excitations for which 9{x, r) is singular at 
the center of the topological excitations. So for this 
type of system, we express the 9{x,t) into two com- 
ponents: 9{x,t) = 9o{x,t) -\- 9i{x,t), where 9q{x,t) is 
the contribution from attractive interaction of the sys- 
tem and 9i{x,t) is the singular part from topological 
excitations. We consider at any arbitrary time r, a topo- 
logical excitation with center at X(t) — {xq[t),tq{t)) 
. The angle measured from the center of topological ex- 
citations between the spatial coordinate and the x-axis 
9i{x^t) — tan^^{ ^" ) . The derivative of the an- 



gle is VxOii 



X{t)) 



, {x - Xo)] 



which has a singularity at the center of the topological 
excitation. Finally we get an interesting result when we 
integrate along an arbitrary curve encircling the topolog- 
ical excitations. J^^dxVx9i{x — X{t)) = 27r. So we 
conclude from our analysis that when a topological exci- 
tation is present in the SQD array, the phase difference 
,0, across the junction of quantum dots jumps by an inte- 
ger multiple of 27r. This topological excitation is nothing 
but the QPS in the (x, r) plane. According to the phase 
voltage Josephson relation , Vj = "jj"^ is the time), 
a voltage drop occurs during this phase slip, which is the 
source of dissipation. This analysis is valid when $ is 
away from ^ . Otherwise system is in the superconduct- 
ing Coulomb blocked phase. 

Now the problem reduces to finding the quantum dis- 
sipation physics of SQD array with modulated Joseph- 
son couplings and critical current. There are a few 
interesting studies, following the prescription of Calde- 
ria and Legget [11|. to uncover the quantum dissipa- 
tion physics of low dimensional tunnel junctions system 
13, lEliSjEI- Our starting quantum action is the same 
with Ref. [13, 13, 14, 1^. We will see the scaling anal- 
ysis of RG equations derive from this action is sufficient 
to explain the experimental findings of Ref. [3-4] . 
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Here, is the standard action for the system with tilled 



wash-board potential [IJ, [13|, [U, |15[ to describe the 



dissipative physics for low dimensional tunnel junction 
systemi, m is the Matsubara frequency. Where Sq = 

la ^(^)^+^(^(^)) ' " = Matsubara frequency 
— ^™ Rq (= GAbkrt) is the quantum resis- 
tance and Rs is the tunnel junction resistance, f3 is the 
inverse temperature. Here V{9) — — |j|cos(^)| -I- 
We would like to exploit the renormalization group (RG) 
calculation for weak potential. Without loss of general- 
ity, we do the analysis for / = case as finite / only 
inclined the potential profile. Since we are interested in 
the low energy excitations, we can ignore the contribu- 
tion of \ujrn\^ compare to \ujrn\- So in Sq, we only consider 
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the second term. 



-/ 7r<I> f'^ 
So = -^\cos( — )\ / cos(9)dT = Vi / cos(d)dT 
2e $0 Jo Jo 

(4) 



We can write the final action as 



- T^El^™!' +^1 f cos{e)dr (5) 

The RG equation of the above action to study the low 
energy excitations is the following: 



dVi 
dlnb 



(1 



1, 



)Vi. 
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We have derived the above RG equation following the 
prescription of Ref. [l6| . Here 5 is a number ratio of 
the two energy scale. The time evolution of the coupling 
constant is Vi{t) = Vi{0)e^^~^^°'^* . We are mainly in- 
terested in the low energy theory of the system, suppose 
we consider the low energy frequency as ujm then the 
corresponding time is i = ln{-^). The coupling con- 
stant at maximum time reduces to Vi = ^i(0)(^)^ " . 
We consider the lowest frequency allowed by the Mat- 
subara allowed frequency quantization, i.e., Um = 2'i:T, 
so Vi(T) (X T~ ~ ^. When we consider a finite sys- 
tem of length, L (i > 1), then we might argue that the 
mode 0{k) is quantized with oji — ^ ^ where uji is the 
lowest frequency and v is the velocity of low energy ex- 
citations. So the effective value of the coupling at the 
lowest frequency is Vi{L) (x L^~~ . We observe that the 
potential Vi increases for a > 1 and decreases for a < 1. 
So a = 1 is the phase boundary. When a > 1 owing 
to the strong dissipation effects particle comes to rest 
at one of the minima of the potential (local or particle 
like character), i.e., the system is confined in one of the 
metastable current carrying state. Hence the system is 
in the superconducting phase. It is interesting to observe 
that dissipation favors to stabilize the superconducting 
phase, a < 1, implies weak dissipation has no effect on 
the potential. The phase fluctuation is large around the 
dots (nonlocal or wave like process). As a result, there is 
no phase coherent state in the system. Therefore there 
is no superconducting phase. 

We also observe from our study that the applied magnetic 
flux has no effect on RG equation at the one-loop level. 
We consider the effect of magnetic flux at the phenomeno- 
logical level, i.e., we replace a by a' = a\cos{^)\. When 
4> is zero or an integer multiple of flux quantum, the flux 
has no effect on the dissipation physics. For the larger 
values of <i>, (small a'), make the quantum fluctuations in 
the system large thereby destroying the phase coherence 
of states. So the higher magnetic field drives the system 
from the superconducting phase to the insulating phase. 
This is consistent with the experimental findings |3| • 
The analytical structure of our derive RG equation (Eq. 



6) has some similarity with the RG equation of single im- 
purity Luttinger liquid [isj . But the initial Hamiltonians 
of these two problems are quite different. 
In the strong potential, tunneling between the minima of 
the potential is very small. In the imaginary time path 
integral formalism tunneling effect can be described in 
terms of instanton physics. We will see that the strong 
coupling physics of our s yste m can be described interms 
of tunneling physics 17 , 13, 3] ■ In the imaginary time 



path integral formalism the potential is inverted and 
therefore the particle can not reside at the maximum of 
the potential for long time and rolls down to one of the 
potential minima. It is convenient to characterize the 
profile of 6 in terms of its time derivative, 

de{T) , 



dr 



Mt 



(7) 



where h{T — Ti) is the time derivative at time r of one in- 
stanton configuration, is the location of the ith instan- 
ton, Ci = 1 and —1 are respectively the topological charge 
of instanton and anti-instanton. Integrated the function 
h from —oo to cxd, dTh{T) = 9{oo) — 9{^oo) = 27r. So 
the appearance of instanton (anti-instanton) is nothing 
but the appearance of QPS in the lumped Squid system 
and it leads to the dissipative phase of the system. Larkin 
et al. [§| have also supported the idea of QPS as a appear- 
ance of instanton (anti-instanton). So our system reduce 
to a neutral system consist of equal number of instanton 
and anti-instanton. Now our prime task is to present the 
partition function of the system. After a few steps, we 
will come to that stage. One can find the expression for 
0{uj), after the fourier transform to the both sides of Eq. 
7 and that yield 



eih{iuj)e^' 



(8) 



Now we substitute this expression for 9(u)) in the second 
term of Eq. 3, finally we get this term as F{Ti — 



Tj)e,ej, where F{Ti - Tj) = ^J2r, 



ln{n " Tj) 



We obtain this expression for very small values of 
( ^ 0). So effectively F{Ti — Tj) is representating 
the Coulomb interaction between the instanton and anti- 
instanton. This term is the main source of dissipation 
physics of the system. Following the standard prescrp- 
tion of imaginary time path integral formalism, we can 
write the partition function of the system 17, 18, 19l |. 



n=0 



Where 

z = e~ 
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the contribution 

{Sinst — V^lC 
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from instanton. 
We may also 



write this expression as 



Z 



n=0 



dT„ 



dTn-l- 



dri 
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(10) 



Here q-r is the auxiliary field which arises during the func- 
tional integral. We again introduce a instead of a to 
consider the effect of applied magnetic field. After ex- 
tensive calculation, we get the partition function, 



Z = 



J Dq{T)e '"'"'^ 



Y,^ I'^lk™!' + 22: drcos(q(T)) 



(11) 

qij) is the auxiliary field. Following the method of pre- 
vious paragraphs , we finally obtain the RG equation 



dz 
dlnb 



{l-a')z 



(12) 



The analytical structure of the quantum action and the 
RG equation is the same with weak potential, it implies 
the following mappings ct' ^ and Vi ^ z. Hence there 
is a duality in this problem between the weak and strong 
potential. One can also find the analytical expression for 
the variation of z as a function of temperature and length 
as we derive for weak coupling case, the only change being 
a' replaced by 1/a'. Fugacity depends on temperature 
and length scale as, 2:(T) (xT"'-^, z{L) cx i^-"'. 
In this complicated system, we estimate the behavior of 
resistance from the behavior of fugacity. It is expected 
to scale as z^, (because the major contribution of volt- 
age/resistance occurs from the second order expansion of 
partition function, i.e., from the square of fugacity). In 
our study, resistance is evolving due to dissipation effect 
at very, low temperature (less than the superconducting 
Coulomb blocked temperature). According to our calcu- 
lations, for large dissipation {a' > 1), R{T) cx RqT^^, 
(3i > 0. Therefore at very low temperature, the system 
shows the superconducting behavior. When a' < 1, the 
resistance of system R{T) oc RqT~^^^, P2 > 0. So at 
very low temperature, the resistance of the system shows 
Kondo-like divergence behavior and the system is in the 
insulating phase. According to our calculations, for large 
dissipation (a' > 1), R{T) oc RqL-^\ 71 > 0. Therefore 
the longer array system shows the less resistive state than 
shorter array in the superconducting phase of the system. 
When a' < 1, the resistance of system R{T) cx RqL^^, 
where 72 > (/3i , /32 , 7i and 72 are independent num- 
bers). So the resistance at the insulating state is larger 
for longer array system than the shorter one. So we find 
the dual behavior of the resistance (voltage) for lower 
and higher values of magnetic field. When a' — 1, i.e., 
$ = (<i)o/7r)cos~^(^), the system has no length scale de- 
pendence superconductor-insulator transition at very low 
temperature. This is the critical behavior of system for 
a specific value of magnetic field. These theoretical find- 
ings are consistent with the experimental observations 
(3^3] ■ Fig. 2 shows the variation of resistance with tem- 
perature. At higher temperature, larger than Tc. SQD 
array system is in the normal phase and the tunneling 



between the dots is the sequential tunneling (one after 
another, tunneling of Cooper pairs) . We have shown in 
Ref. [201, that superconductivity occurs due to the co- 
tunneling effect (higher order tunneling of Cooper pairs, 
virtual process). A presence of finite resistance at the 
superconducting phase (between Ti and T2 ) due to the 
dissipation effect and also for the presence of finite tun- 
neling conductance (i.e, the finite resistance). Low resis- 
tance superconducting phase or insulating phase occurs 
at very low temperature, smaller than the superconduct- 
ing Coulomb blocked temperature. 

Low resistance superconducting phase or insulating 
phase occurs at very low temperature (it will occur 
for few mili-kelvin) , smaller than the superconducting 
Coulomb blocked temperature. In experiment, they have 
measured up to 50 mK, so they have not found the decay- 
ing tendency of resistance at very low temperature. We 
have not considered the the classical phase {Ej >> Ec) 
of the system. In this phase, one can also obtain the dis- 
sipative phase with phase slip centers, but one can loses 
the informations of intermediate quan tum phases for the 
system Chakarvarty et al. 13 1 and Larkin et al. 

[9| had studied QPS for the classical phase. Currently 



Fistual et al. [21| have done some interesting work on 



collective Cooper pair transport in the insulating state 
of one and two dimensional Josephson junctions array. 
They have studied the current-voltage characteristics re- 
vealing thermally activated conductivity at small volt- 
ages and threshold voltage depinning. Our analytical 
approach to study the one dimensional mesoscopic squid 
system is quite different from them [ilj. 
We have studied the clean system. Here we explain the 
effect of impurities in superconducting dots or in the tun- 
nel barrier. Nonmagnetic impurities were found not to 
affect the Josephson supercurrent. This is in agreement 
with Anderson's theorem [2^ of dirty superconductor. In 
presence of paramagnetic impurities, there is an exchange 
interaction between the spin of conduction electrons with 
the magnetic impurity spin. The spin of the magnetic im- 
purity polarize the spin of electrons and interferes with 
their tendency for pair formation in the singlet state. As 
one increase the probablity (F) of scattering with spin 
flip. For highervalues of F, the system enters into the 
gapless region and above a critical value ( Fc = , 
A(0, 0) being the order parameter for the superconductor 
at zero temperature in absence of impurity. ) the super- 
conductivity is destroyed (23| . For such situation, our 
system is in the Coulomb blocked insulating phase. The 
Josephson supercurrent is also zero for a clean dot when 
the applied magnetic flux of the system is half-integer 
multiple of flux quantum. At those values of magnetic 
flux, our system is in the Coulomb blocked insulating 
phase. 
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3. CONCLUSIONS 

We have observed the length scale dependent 
superconductor-insulator quantum phase transition and 
the dual behaviour of the system for smaller and larger 
values of magnetic field, in a one dimensional mesoscopic 
lumped Squid systems. We have also observed a critical 
behaviour where the resistance is independent of length 
at a critical magnetic field. We find that weak and strong 
potential results are self dual. Our theoretical findings 
have experimental relevance of lumped Squids system 
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